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I. BACKGROUND 



A. INTRODUCTION 

According to Bartholomew and Forbes [Ref. 1: p. 1], 

Manpower planning is often defined as the attempt to 
match the supply of people with the jobs available for 
them. 

This is a concept that is an integral part of military force 
planning, today. All branches of the Department of Defense 
are involved with studies concerning manpower planning. 
Since large numbers of personnel are involved, the use of 
manpower models to determine the results of changes to the 
existing system before they are actually implemented can 
prevent costly future problems. This is especially true in 
terms of having the right number of personnel in the rignt 
places with the right grade levels and skills, at the right 
time. To quote again from Bartholomew and Forbes [Ref. 1: 

p. 61 , 

A manpower model is a mathematical description of how 
change takes place in the system. 

Much work has been done in the area of manpower/ 
organizational modeling. The text by Bartholomew and Forbes 
(referenced above and hereafter referred to as BaF) compiles 
information cn the subject from numerous disciplines and 
standardizes notation to describe some of the basic models 
and statistical techniques used in manpower planning. A 
central theme of this work is the idea of depicting an 
organization as a system of sto cks and flows [Ref. 1: pp. 

3,4]. The personnel within the organization under study are 
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divided into classes, or cat 
istics (for example: age. 

The numbers of personnel w 
time are considered stocks, 
personnel from one catego 
interval of time, such as a 
Various mathematical mo 
from one year (or other t 
described in EaF. Of parti 
sition model that is based 
In addition to describing th 
provides a BASIC computer pr 
calculate future stock vecto 
into an APL program at the 
and is available on the N 
system 370. The APL program 
stocks from cne year to th 
initial stocks, recruitment 
and the transition probabili 
be done repetitively and th 
stocks forever. 

There are two versions 
that are well developed in 
flow of recruits into the s 
total system size. Also of 
options) that deal with a 
either recruitment or total 
or multiplicative fashion, 
computation of future stocks 

(1) recruitment is fix 
pi icat ively; 

(2) the total system 
plicat ively. 



egories , by specific character- 
time in service, grade, etc), 
ithin categories at a specific 
Flows are the movements of 
ry to another during a unit 
year. 

dels that predict future stocks 
ime period) to the next are 
cular interest here is the tran- 
on the theory of Markov Chains, 
e model [Bef. 1 ; pp. 85-98], BaF 
ogram [Bef. 1: pp. 248-260 ], to 
rs. That program was converted 
Naval Postgraduate School (NPS) 
PS mainframe computer, an IBM 
, like its predecessor, computes 
e next based on input such as, 
and wastage (or attrition) rates 
ties for the stocks. This can 
eoretically could compute future 



of the basic Harkov Chain model 
BaF. One is based on a fixed 
ystem and the other on a fixed 
interest are other versions (or 
system undergoing growth where 
system size grows in an additive 
The APL program allows the 
for the cases where: 
ed, grows additively or multi- 
size grows additively or muiti- 
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Steady state is the condition which a system reaches, 
after sufficient time has passed, so that the stocks no 
longer change- There have been steady state models devel- 
oped for the cases where there is a fixed flow of recruits 
into the system and where the total system size is fixed. 
In both these cases the number of personnel in each category 
remains fixed after steady state is reached by the system. 
The vector whose components are these numbers is called the 
steady state stock vector (S.S.S.V. ). 

Under conditions of growth, the steady state of a system 
must be redefined. In those cases, steady state does not 
exist in the above sense, because the system is undergoing 
constant expansion and so stocks will never reach a fixed 

size. However, it turns out that the stock sizes do reach 

constant proportions and so it is possible to talk about a 
steady state distribution vector (S.S.D-V-). Components of 
this vector represent the pro po rtions of personnel that 
remain in each category forever once steady state has been 
reached. In effect, in taose cases where the distribution 
of personnel among the categories no longer changes with 

time the system has reached a steady state "in 

dis tr ib ution. " 

Although, at the present time, the APL program at the 
Naval Postgraduate School predicts future stocks for a 
system that grows under conditions mentioned above, it is 

able to supply a steady state stock vector (S.S.S.V.) only 
for the case where there is a fixed flo w of recru its into 
the system. The only way to get the steady state distribu- 
tion vector under the various conditions of growth is to 
continue to run the program, i.e., to forecast stocks for as 
many years into the future as necessary until steady state 
is reached. This is not a satisfactory method for several 
reasons : 
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(1) it takes too long to run the model when the system 
reaches steady state only after many years have 
passed , 

(2) it is a trial and error process of going back and 
forth, in time, to try and find out at what point 
in time steady state has been reached, 

(3) occasionally, even if the program is run out as 
far as the present formatting of the pregram 
allows, there can still be instances where steady 
state has not yet been reached. 

This thesis uses the ’basic prediction equation’, as 
described in BaF [ Bef . 1: pp. 86-88 ], to correct the afore- 
mentioned deficiencies by developing models to compute 
analytically the steady state distributions for systems 
under varying conditions of growth. Also developed was an 
APL program, to oe added as a subroutine to the main 
program, that computes the S.S.D.V. for these options. 

The next section in this chapter sets out the necessary 
notation and equations to enable the reader to follow the 
development of these steady state models. For more amplifi- 
cation on the casic Markov model see BaF [Ref. 1: pp 3-8, 
86-132]. 

B. NOTATION AND BASIC EQUATIONS 

Stocks, representing personnel at an organization, are 
divided into k categories. The stock vector 

2 (t) = ( n, (t) ,r. 2 (t) , . . . ,n K <t) ) 

is a row vector where each element n^ (t) , represents the 
number of personnel in category i at time t. Time is meas- 
ured discretely with a time interval of unit length denoted 
(t-1,t) . The most frequently used period of time is the 
year although it is sometimes convenient, or necessary, to 
use quarters, months, etc. 
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Flows describe the movements of personnel. There are 
two basic types of flows in a system. The first type of 
flow is internal to the system and can be thought of as 
transf ers from one category to anoth'er within the system. 
Demotions and promotions would be examples of this type of 
flow. These are denoted by the k by k matrix, N(t-1), where 
each element, n ( *j (t-1), represents the number of personnel 
moving from category i to category j during the time 
interval (t-1,t). 

The other type of flow involves the two-way transfer 
between the system (organization) and the outside world. 
First, there are w as ta ge (or attri t ion ) flows of the number 
of personnel leaving the system from category i, during the 
time interval (t-1,t), denoted by n^ K+ .(t-1), (i = 1,2, ...,k 
and t = 1,2,...). Then there are r ecr uit ment (or ac cess ion) 
flows of personnel into category j from the outside world, 
during the period (t-1,t) , denoted by n G j (t) , (j = 1,2,...,k 
and t = 1,2,...). Following BaF [Ref. 1: pp. 3,4 ], the 

notation for recruitment flows is slightly different from 
that for internal flows and wastage flows. This is to 

remind the reader that in most organizations attrition and 
internal flows (transfers) are accounted for first and then 
recruitment is determined at the end of the time period. 

The Markov Chain model uses flow rates (or probabili- 
ties) rather than actual numbers. P, a k by k matrix of 
probabilities has elements 



p. - . t i,j = 1,2,... ,k . 

n L (t-1 ) 



There are certain basic assumptions 
model needs tc meet to be valid: 

(1) the transition probabilities, 
over time; 



that a Markov Chain 




do not change 
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(2) the population in each category is homogeneous so 
that p^ represents the probability of each indi- 
vidual in category i moving, independently of any 
other individual, to category j. 

While these assumptions may not always be met, there 
have been enough studies done using the model to show they 
are sufficiently realistic to justify the use of the model. 
Even in cases where these assumptions don't hold very well, 
the model was fcund to be quite useful for predicting future 
stocks [Refs. 2,3]. 

Next, w = (w., ,w 2 ,...,w K ) is defined as the wastage 

(attrition) rate vector where each element, w- , is defined 
as 

u . = 

L n- (t- 1) 

Ncte that, 

K 

Pj . ^ ^ — 1 , i - 1,2,««.,k . 

jsl ) 

Finally, r = (r 1 ,r 2 , . . . , r K ) is defined as the recruit- 
ment proportion vector where each element 

r . = 

J E (t) 

is the proportion of total recruitment that will be allo- 
cated to category j, with 

K 

R <t) = J n oj 

j = i J 

being the total recruitment during the time interval 

( t— 1 , t ) . This implies ^ r- =1. 

j 1 ' 

BaF [Ref. 1: pp. 87,88], snow that future stocks may be 
calculated using the equation 
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♦ fi ( t) r- 



t 



(eg a 1.1) 



n- (t) = 



I 



L- I 



n t (t-1) 



C J 



for j= 1,2, ...,k and t = 1,2,... . This equation states, 

that the stocks in category j at time t are equal to the 
number of personnel that moved into category j from anyv?here 
within the system during time period (t-1,t), plus the 
number of personnel recruited into category j during the 
same period. In matrix form the equation is 



n (t) = n ( t- 1) P + E (t) r t = 1,2,... 



(ecn 1.2) 



BaF [Hef. 1,: p. 88] refers to this as "the basic prediction 
equation." 

Another equation dealing with changes in total system 
size is also developed in BaF [Bef. 1: pp. 94-96 ]. To see 
this, note that 

K 

N(t) = 2n c (t) 

L- l 

is the total system size at time t and 

M(t) = N ( t ) - N (t- 1 ) (eqn 1.3) 



represents the change in the total system size during the 
time interval (t-1,t). This implies that total recruitment 
during (t-1,t) can be expressed as: 



K 

B (t) = ^n ( ( t— 1 ) w. + M (t) > 0 

i 3 f 



(eqn 1.4) 



The first term in this expression is recruitment that is 
done to replace those who left (attrition from the various 
categories) while the second term consists of recruitment 
that is due to the change in total system size. 

If M (t) =0, then only those personnel who leave the 

system are reolaced, which implies that the total system 
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size remains fixed. If M (t) > 0, new personnel are supplied 

to cover an increase in total system size, i.e., to fill new 
jobs. If M { t) < 0, then some jots must be eliminated as 
soon as they are vacated. Of course, E (t) i 0 always. 

Referring back to Equation 1.1, replace R (t) using 
Eguation 1. 4 and let 



g 



* « 

‘J 




+ w • r • . 

6 J 



Then, as Ba F [Ref. 1: pp. 94,95], 
K 

n i M = y n L Ct-1 ) q.- + M (t)r . 

J £ C cj J 

for j= 1,2, ...,k and t = 1,2,... 

notation let 



shows 



, (egn 1.5) 

. To put this into matrix 



Q — [ q ' • } r i»j — 1,2,...,k , 
‘•J 

or 



Q = P + w 'r . 

Then Eguation 1.5 may be expressed as 

n_(t) = _n(t-1)Q + M(t)£ , t = 1,2,... (egn 1.6) 

Equations 1.2 and 1.6 form the basis from which five 
steady state submodels (or options) will be derived. The 

next section introduces these options and the direction of 
develop ment . 

C. OBJECTIVES 

The objective of this thesis is tc develop steady state 
equations modeling five different- conditions of growth. 
Each of these conditions will be addressed separately as a 
"growth opticn." The first three options deal with 
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different types of recruitment. They use Equation 1.2 where 
R (t) is the total recruitment during the time period (t-1,t) 
and in these three options, R (t) is the only term changing. 
The last two options deal with changes in total system size. 
They utilize Equation 1.6 where. M(t) represents the change 
in total system size during the period (t-1,t) and in these 
two options, M (t) is the only term that changes. A complete 
list of these five options is: 

(1) Fixed recruitment: R (t) = F, t>1 , R>0 fixed, 

i.e., recruitment is always constant; 

(2) Additive increase (decrease) of recruitment: 

R (t) = R + (t-1) M , t>1 

where M is the constant amount of yearly increase 
(decrease) ; 

(3) Multiplicative increase (decrease) of recruitment: 
E(t) = Re*-i, t>1, 

where 100(9-1) is the percent of yearly increase 
(decrease) ; 

(4) Additive increase (decrease) in system size: 

M ( t) = M , t> 1 

where M is the constant amount of yearly increase 
(decrease) in total system size; 

(5) Multiplicative increase (decrease) in system size: 

M ( t ) = (0-1 ) 0*~ 1 N (0 ) , t>1 

where 100(9-1) is the percent of yearly increase 
(decrease) in total system size. 

If H =3 in option 2 or 6 = 1 in option 3, then those 
options revert to option 1, of fixed recruitment. If M = 0 
in option 4 cr 6 = 1 in option 5, then M(t) =0 and the 
result is a system with a fixed total system size. That 
case will be dealt with as subcases of options 4 and 5. The 
fixed total system size case and option one, where R is 
fixed, are the only cases where it is possible to achieve a 
steady state stock vector, i.e., where the stocks don't 
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change with time. In all other cases steady state will be 
considered to have been reached when a steady state distri- 
bution of personnel is achieved. 

Chapter II addresses, in depth, the analytical develop- 
ment of the steady state equation for each of the options. 
It starts with an account of the work done on option 1 in 
BaF. Then it proceeds with the derivation of the steady 
state equations for options 2 through 5 from either Equation 
1.2 or Equation 1.6 . It also presents examples, one a 
four-grade hierarchical system, the other a three grade 
non-hierarchical system. In addition to formulating analyt- 
ical models fcr options 2 through 5, an APL program has also 
been developed that will compute the steady state distribu- 
tion tor those options. This program is contained in 
Appendix F. Appendices A through E are computer printouts 
of terminal sessions of the two examples for each of options 
1 through 5. 
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II. OPTION DEV EL OPMENT 



The options mentioned in Chapter I will be discussed 
here under the usual restrictions of manpower systems. 
These restrictions imply that w; > 0 for all i=1,2,...,k, 

since in manpower systems attrition must be allowed from any 
category- This means that the P matrix, as defined in 
Chapter I, is composed of non-negative elements and has row 
sums all strictly less than one. The Q matrix, also defined 
in Chapter I, is a stochastic matrix (i.e. , it is composed 
of non- negative elements and has row sums all equal to one). 
These facts along with the following two mathematical theo- 
rems are needed in the derivation of the steady state 
results for the five options. 

Theorem _1 

If P is a matrix composed of non-negative elements and 
row sums strictly less than one and 0 is a scalar then 
the matrix 0I-P has a unique inverse for all values of 
e > 1 [Bef. 4: p. 42 ]. 

Theore m 2 

If Q is a stochastic matrix and 6 is a scalar then the 
matrix 0I-Q has a unique inverse for all values of 
9 > 1 [ Ref. 4: pp. 60,6 1 ]. 



A„ RECRUITMENT OPTIONS 

The following three options are all variations of 
Equation 1.2, where R (t) , representing total recruitment, is 
the term that takes on different forms. Consequently, the 
steady state equations for these options are also obtained 
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in somewhat different forms. These steady state equations 
are derived in the next three subsections. 

1. Opt ion _1 : Fixed Recruitmen t 

The steady state equation for this first option is 
developed in BaF [Ref. 1: p.90]. The analysis presented 

here is a brief summary of that in BaF, using the same 
assumptions. Recalling the basic prediction equation, 
(Equation 1.2), and assuming that 

limit R (t) = R is fixed 

t — oo 

it follows that the steady state stock vector, 

n = limit n (t) 
t — oo 

also exists [Ref. 4: pp. 40-43,48-50 ]. Therefore, letting 

t— oo in Equation 1.2 results in 

n = nP + Rr , 



or 



_n (I - P) = Rr . 

Finally, knowing the matrix I-p has a unique inverse (see 
Theorem 1), the result 

n = Rr <1 - P) ~ l , (eqn 2. 1) 

is obtained for tne steady state stock vector (S.S.S.V. ). 
After steady state has been reached, it follows that 

£(t) = n_ , 

for "large” t. The following two examples will be used to 
illustrate the use of this steady state equation. 
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Example (a) 

This first example uses data from Example 4.6 in BaF 
[ Bef . 1: p. 97], which is taken from the women officer's 



system of one of the British services, 
ical, four category system. 

The inital stock vector is: 

n_(0) = .(129, 74, 28, 1 1) 

The transition matrix is: 

.102 

.83 
0 
0 



It is a hierarch- 



P = 





The total recruitment is: 

B(t) = 35 , fixed for all t>1 
The recruitment proportion vector is: 

r = (1, 0, 0, 0) . 

Then 



I - P = 





and the inverse is 



(I - P)~i = 



676 


2.206 


.763 


. 257 


0 


5. 882 


2. 034 


.685 


0 


0 


7. 519 


2.532 


0 


0 


0 


10. 204 
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Using Eguation 2.1 



n = (35) (1, 0, 0, 0) 




2.206 


.763 


. 257 


5.882 


2.034 


. 685 


0 


7.519 


2. 532 


0 


0 


10. 204 



= (128. 66, 77. 21 , 26. 705, 8.995) . 



after rounding to the nearest integer the S.S.S.V. is: 
n = (129, 77, 27, 9) 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 
can be found in Appendix A, confirming the steady state 
stock vector above. As the printout shows, steady state is 
achieved by t = 10. 



Example (b) 

This is a n on-hierarchical, three-grade system. 
The inital stock vector is: 



H(O) = 


(300, 


200, 


100) . 


matrix 


is: 








/ . 6 0 


. 15 


.05 \ 


P = | 


• 75 


. 20 


0 




\.05 


0 


.90 / . 



The total recruitment is: 

R(t) = 100 , fixed for all t> 1 
The recruitment proportion vector is: 

r = (.70, .07 , . 23) . 
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Then 



I - P = 



and the inverse is 



(I - P)-i = 



. 4 


-. 15 


-.05 \ 


.75 


. 8 


0 


. 05 


0 


,10 / 


4.267 


.8 


2. 133 


4 


2 


2 


2. 133 


. 4 


11.066 



Using Equation 2.1 



4. 267 .8 

n_ = (100) (. 70, .07, .23)1 4 2 



2. 133 . 4 1 1.066 




= (375.7, 75. 2, 417.8) . 

After rounding to the nearest integer the S.S.S.V. is: 
n, = (376, 79, 4 18) 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 

can be found in Appendix A , confirming the steady state 
stock vector above- As the printout shows, steady state is 
achieved by t = 50. 



2. Opti on 2: A dditive Re cru i tmen t 

This option addresses the fixed additive increase 
(decrease) in recruitment for each time period. Here, 
R + (t — 1) n , replaces R (t) in Eguation 1.2 resulting in 



n.( t) = n (t- 1) P + Rr + ( t— 1 ) Mr. 



(eqn 2.2) 
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In this option, M > 0 must be true, otherwise there 
would be negative recruitment (R (t) <0) which does not make 
sense in manpower models. This means there are two possible 
cases. 

Case 1: M = 0 

Then 

limit R (t) = R , 

t — OO 

which is identical to option 1. 

Case 2: M > 0 

Assume that, 

limit = n_ 

t — tM 

exists. This is the same as assuming that for "large" t, 
ri(t) behaves like tMrT . Using the above assumption. 
Equation 2.2 can be rewritten as 





n (t) 




n(t-1) 


n 


R (t- 1) M 






tM 




tM 


r 


▼ IT • 

t M tM 


* JC 


or 


n (t) 




n<t-1) 




(t- 1 ) _ R 


(t-1) M 

i 




tM 




(t-1) M 




• P + r 

t t M 


tM 


Tak ing 


the liiti 


t as t— » 




results in 








A 

n 


II 

\*> 




r 




or 




A 


(I-P) = 


£ 


- 




Since 


the I-P 


matrix has 


a 


unique. inverse (se 


e Theorem 



the result 

n* = r (i-P)-i 
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is obtained. The steady state distribution vector is 
achieved by normalizing n|, which results in 



yN 

n 



n = 



nl* 



L [I 

nl • 



(eqn 2.3) 



After steady state has been reached it follows that 
the formula 

n(t) = tun 1 for "large" t, 

can be used to approximate the stocks at time t. The two 
examples from section Al are worked out below to illustrate 
Equation 2.3 . 

Example (a) 

This is a hierarchical, four category system. 

The inital stccx vector is: 

_n(0) = (1 29, 74, 28, 1 1) 

The transition matrix is: 



P = 




, 102 
,83 
0 
0 



0 

,0 46 
, 867 
0 




The total recruitment is: 

R (t) = P. + (t- 1 ) M where R = 35 and M = 5 - 

The recruitment proportion vector is: 

r = (1, 0, 0 , 0) 
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Then 



I - P = 





and the inverse is 



(I - P)~i = 





Using Equation 2.3 



n = 



3.676 2.206 



,763 



.257 



% (1» Or 0, 0)/ 


0 


5.882 


2.034 


. 685 \ 


Itl* 6.902 1 


0 


0 


7.519 


2. 532 1 


V 


0 


0 


0 


10.204 / 



= (. 53, . 32 , . 11 , .04) 



is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, 100, 

900, and 975 can be found in Appendix B, confirming that in 
steady state there will be 5 3 5?, 325?, 11%, and 4% respec- 
tively, in categories 1,2,3, and 4. As the printout shows, 

steady state distribution (to two decimal accuracy) is 
achieved oy t = 975 . 



Exa mple (b) 

This is a non-hierarchical, three-grade system. 
The inital stock vector is: 

n (0) = (300, 200, 100) 
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The transition matrix is: 



60 .15 .05 

P = | .75 .20 0 

05 0 .90 

The total recruitment is: 

R (t) = R + (t- 1 ) M where R = 100 and M = 25 

The recruitment proportion vector is: 

r = (.70, .07, .23) . 



Then 



I - P = 




-. 15 -.05 



8 

0 




and the inverse is 



(I - P) 








Using Equation 2.3 









/ U - 


267 


8 


2. 


133 




(.7, 


.07, . 


23) / 


0 2 


.0 




2.0 


n - ~ 






4. 




“ ~ n1_’ 


8. 


728 


1 
















\ 2. 


133 


4 


1 1 . 


066 


= ( • 


43 , .09 


i, .48) 












is the 


stea d y 


state 


distr ibu 


ticn ve 


cto r 


m 






The computer 


printout of 


the 


ter 


min 


this ex 


ample , 


showing the 


stocks 


for 


yea 


rs 



50, 



for 

100 , 
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350, and 400 can be found in Appendix B, confirming that in 
steady state there will be 43%, 9%, and 43% respectively, in 
categories 1,2, and 3. As the printout shows, steady state 
distribution (to two decimal accuracy) is achieved by 
t = 400 . 



3- O pt ion 3: Multi p l ic ativ e Be er uit ment 

This cption addresses the multiplicative increase 
(decrease) in recruitment for each time period. Here, 
Re*- 1 , replaces E (t) in Equation 1.2 resulting in 



n(t) = n (t- 1) P + E0 t -»r 



(eqn 2.4) 



The values that 6 can assume play an important part in the 
derivation of the steady state equation. First, © >0 must 
be true, otherwise in even numbered years there would be 
negative recruitment ( R(t)<0 ) which does not make sense in 
manpower models. There are three cases to be examined. 

Case 1 : 0< 8 <1 



Then 



limit R (t) = 0 , 

t — ► oo 



because e*- 1 — ►O when 0< € <1 . Therefore, letting 

t— *> in Equation 1.2 results in 



or 



n = nP , 



n ( I - P) = 0_- . 

Since I-P has a unique inverse matrix (see Theorem 1) , the 
only solution to the above equation is n = 0.. This means 

that such a system becomes empty in steady state. 
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Case 2: 0 = 1 

Then 

limit fi (t) = R , 

t — •* oo 

which is identical to option 1. 

C as e 3 : 0 >1 

Assume that, 

limit — = n 

t — oo e* — 

exists. This amounts to assuming that for "large" t, n{t) 
behaves like e* - £ . Using the above assumption. Equation 
2.4 can be rewritten as 

n (t) iL(t-l) Re*-* 

=— = — *P + 1 r 

e r e*-i 0 0 

Taking the limit as t — ► results in 

^ ^ 1 n B 

n = n« - *P ♦ - r 

“ 0 9 ”, 



or 

[ei-P] = Hr . 

Finally, knowing that the matrix ©I-P has a unique inverse, 
since 0 > 1 (see Theorem 1) , the result 

n = Er [ SI— P ]-* 

is obtained. When n is normalized, the resulting vector is 
n, the steady state distribution vector where. 



_ RrJ ei-P ]- 1 (egn 2.5) 
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After steady state has been reached it follows that 
the formula 



B(t) 



■fc A 

0 n 



for "large" t. 



can be used to approximate the stocks at time t. Two exam- 
ples are worked out below to illustrate Equation 2.5 



E xample (a) 

This is a hierarchical# four category system. 
The inital stccK vector is: 

n (0) = (129# 74, 28, 1 1) 

The transition matrix is.: 




The total recruitment is: 



£ <t ) = E0^-i where R = 35 , 6 = 1.0 1 . 

The recruitment proportion vector is: 



jc — ( 1 , 0 , 0 , 0 ) . 

Then 




and the inverse is 



28 



[ ei - p]-i 



Using Equation 2.5 



546 


2.009 


.646 


.198 


0 


5. 556 


1.787 


.546 


0 


0 


6.993 


2.137 


0 


0 


0 


9. 259 



n = 



A 

n_ 

nl ' 



(35) ( 1, Q f 0, 0) 
223.965 



/ 3. 54 6 


2. 009 


.646 


.198 \ 


f 0 


5. 55o 


1.787 


.546 1 


0 


0 


6.993 


2.137 1 


\ o 


0 


0 


9.259 / 



= (- 554, .3 14, . 101, . 031) 

is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 

can be found in Appendix C, confirming that in steady state 
there will be 55%, 31%, 10%, and 3 % respectively, in catego- 
ries 1,2,3, and 4. As tne printout shows, steady state 
distribution (to two decimal accuracy) is achieved by 
t = 50 . 



Example (b) 

This is a n on-hierarchica 1, three-grade system. 
The inital stock vector is: 

n_(0) = (300, 200, 100) 

The transition matrix is: 

( .60 -15 .05 

. 75 . 20 0 

.05 0 .90 
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The total recruitment is: 



R (t) = Re <_l where R = 100 
The recruitment proportion vector is: 

r = (.70, .07, . 23) . 

Then 



© = 1.01 . 



61 - P 



and the inverse is 



[01 - p]-i 




( 4 . 03 
3.73 

1.33 



-.15 
. 8 1 

0 



.75 

1.92 

. 34 




1. 83 
1.69 

9. 92 




Using Eguaticn 2.5 



n = 





/ 4. 0 3 


.75 


1.83 


n (100) (.70, .07, .23) | 


3.73 


1.92 


1.69 


Hi* " 792.15 


\ 1.83 


.34 


9.9 2 



= (-44, . 09, . 47) 

is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 

can be found in Appendix C, confirming that in steady state 
there will be 44%, 9%, and 47% respectively, in categories 
1,2, and 3. As the printout shows, steady state distribu- 
tion (to two decimal accuracy) is achieved by t = 100 . 
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B. TOTAL SISTEH SIZE OPTIONS 



The options in this section are variations of Equation 
1.6 where !1 (t) , representing the change in total system 
size, is the term that will take on different forms. As a 
result, the steady state equations for these options are 
also obtained in somewhat different forms. In the next two 
subsections the steady state equations will be derived for 
the system size options. 

1- O ptio n 4 : Additive Syste m Size 

This option addresses the additive increase 
(decrease) in total system size for each time period. In 
this case, a (t) = M, a fixed amount of increase (decrease) 

per time period. Rewriting Equation 1.6 results in 

_n(t) = n_(t- 1) Q + Mr_ . (egn 2.6) 

There are 3 cases to be considered. 

Case 1: M < 0 

In this case jobs are being eliminated by a fixed 
amount, M, each year which implies that eventually N (t ) = 0 
is reached in a finite amount of time. This implies that 
n = 0_ is the steady state stock vector. This means that 
sucn a system becomes empty in steady state. 

Case 2: M = 0 

When M(t) = 0, N (t) is fixed and can be denoted N. 
It follows that Equation 1.6 can be written as 

(eqn 2 . 7) 



n (t) = n_(t- 1) Q . 

Since the 



limit n(t) = n 

t— oo 



31 



exists [Bef. 4: pp. 40-43, 48-50 ], Equation 2.7 can be 

written as 



Because the first column of I-Q is a linear combination of 
all the other columns the inverse of I-Q does not exist. 
So, to solve for n, it is necessary to introduce the addi- 
tional constraint of n1_' = N [Ref. 5: p. 80 ]. Replacing 

the first column of the I-Q matrix with a column of ones, 
denoted by (I-Q)*, and replacing the first element of JO with 
N, denoted _0 * , results in the following equation 

_n (I-Q) * = 0* 

Taking the inverse results in 

_n = 0*[ (I-Q)*]” 1 , (egn 2.8) 

the steady state stock vector. After steady state is 
reached, then £{t) = _n for "large" t. Two examples are 

worked cut below to illustrate Equation 2.8 . 

Exampl e (a) 

Thxs is a hierarchical, four category system. 

The inital stccx vector is: 



n = n Q , 



or 



n (I-Q) = 0 



n(0) = (129, 74, 28 , 1 1) 



The transition matrix is: 



P 



( 



728 .102-0 0 

0 .33 .046 0 

0 0 .867 .033 

000 .902 
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and 



Q = 



898 


. 102 


0 


0 


124 


.83 


.0 46 


0 


1 


0 


.8 67 


.033 


098 


0 


0 


.902 



The 0* vector is: 

0* = ( 242, 0, 0, 0 ) 
The recruitment proportion vector is: 

r = (1, 0, 0, 0) 

Then 



I - Q = 



/ . 102 


-. 102 


0 


0 \ 


/ -.124 


. 17 


-.046 


0 l 


I-' 


0 


. 133 


-.033 I 


\ 098 


0 


0 


. 098 / 



first column 


w ith 


ones 


g ives 




/ 




•. 102 


0 


0 \ 


II - Q)* = I 


1 


. 17 


-. 046 


° I 


1 


0 


. 133 


-.033 / 


\ 


1 


0 


0 


. 098 / 



and the inverse is 



((I - Q ) *)“ 1 = 



.533 


.320 


- 1 1 1 


.037 


-4.582 


3.133 


1.084 


. 365 


-5.354 


-3.212 


- 6.408 


2.158 


-5.435 


-3.261 


-1.128 


9. 824 



Using Equation 2.8 
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n = (24 2, 0, 0, 0) 



r .533 


.320 


.111 


.037 


-4.582 


3. 133 


1.084 


. 365 


-5.354 


-3.212 


6.408 


2. 158 


-5.435 


-3.261 


-1 .128 


9. 824 



= (128.99, 77.44, 26. 86 , 3.95) . 

After rounding to the nearest integer the S.S.S.V. is: 



n_ = (129, 77, 27, 9) . 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 

can be found in Appendix D, confirming the steady state 
stock vector above. As the printout shows, steady state is 
achieved by t = 50 



Example (b) 

This is a n on-hierarchical ,' three-grade system. 
The inital stock vector is: 

n (0) = (300, 200, 100) 

The transition matrix is: 

( .60 .15 .05 

.75 .20 0 

.05 0 .90 

and 

( . 74 . 164 .096 

.785 .2035 .0115 

.085 .0035 .9115 

The 0* vector is: 

0* = ( 600, 0, 0 ) . 
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The recruitment proportion vector is: 



Then 



r = (.70, .07, .23) . 



/ . 26 


-. 164 


-.096 \ 


I - Q = (-.785 


.797 


. --012 ] 


\-. 085 


-. 004 


. 089 / 



Replacing the first column with ones gives 





n 


-. 164 


-.056 


(I - Q) * = 


h 


.797 


-.012 




\ i 


-.004 


.089 



and the inverse is 



((I - Q )*)" 1 



/ .430 


.090 


.479 \ 


1 -.611 


1.127 


-.516 ) 


\-4.888 


-.581 


5.869 ' 



Using Equation 2.8 



( . 430 
-.611 

-4.888 



.090 .479 
1.127 -.516 
-.981 5.869 



= ( 258.0, 54.4, 287.4) . 



After rounding to the nearest integer the S.S.S.V 
n = (258, 54, 287) . 



. is 
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The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 ana 100, 

can be found in Appendix D, confirming tne steady state 
stock vector acove. As the printout shows, steady state is 
achieved by t = 50 . 

Case 3: M > 0 

Equation 1.6 can then be rewritten as 

n (t) n (t- 1) M 

= Q + r 

N (t) N (t) N (t) “ . 

Iterating the eguation, N(t) = N <t — 1 ) + M, for t = 1,2,... 

will shew that, N (t) = N (0) + tM for all t = 1,2,... and 

the refo re, 

M N ( t- 1 ) 

limit =0 and limit = 1 

t — oo N(t) t— <*» N(t-1) + K 

Then writing 



n (t) n_ < t- 1 ) N (t- 1 ) H 

= • • Q + r 

N (t) N (t— 1 ) [N(t-1) + a] N(t)“ , 



and taking the limit as t — ► oo yields 



n = b Q , or _n (I-Q) = 0_ 

where 

. n(t) 

n = limit 

t— oo N (t) 

is the steady state distribution vector. 

This is identical to the M = 0 case, except here the 
constraint n1_* = 1 must hold. Therefore, by replacing the 
first column of the I-Q matrix with a column of ones. 



36 



denoted by (I-Q)*, and replacing the first element of 0 with 
1, denoted _0**, the solution is 

n_ = £**[ (I-Q) * , (eqn 2.9) 

the steady state distribution vector. 

After steady state has been achieved, _n(t) can be 
found using 



n_(t) = Q[ N (0) + ti1]n 

for "large" t. Two examples are worked out below to illus- 
trate Eguation 2.9 . 



Example (a) 

This is a hierarchical, four category system. 
The inital stock vector is: 



il(0) = (129, 74, 28, 1 1) 
The transition matrix is: 



P = 





and 



Q = 




102 

,83 

0 

0 



0 

.046 
.8 67 
0 




The 0** vector is: 



0 ** = ( 1 , 0 , 0 , 0 ) 
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The recruitment proportion vector is: 



r = (1, 0, 0, 0) 



Then 





/ . 102 


-. 102 


0 


° \ 


I - 


/ -.124 

2 = 


. 17 


-.046 


0 1 




r * 1 


0 


.133 


-.033 / 




\-. 098 


0 


0 


. 098 / 


Replacing the first 


column with 


ones , 


gives 





(I - Q)* = 




and the inverse is 



((I - Q ) *)~ 1 =| 



.53 3 


.320 


.111 


. 037 


-4.582 


3. 133 


1 .084 


. 365 


-5.354 


-3.212 


6.408 


2. 1 58 


-5.435 


-3.261 


-1.128 


9. 824 



Using Eguation 2.9 



n = (1, 0, 0,0) 



.533 


-320 


.111 


. 037 


-4.582 


3-133 


1.084 


. 365 


-5.354 


-3.212 


6.408 


2. 158 


-5.435 


-3-261 


-1 .128 


9. 824 



= (.53, . 32, . 11 , . 04) 



is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, 100, 

900, and 950 can be found in Appendix D, confirming that in 
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steady state there will be 53%, 32%, 11%, and 4% respec- 
tively, in categories 1,2,3, and 4. As the printout shows, 

steady state distribution (to two decimal accuracy) is 
achieved by t = 950 

Example (b) 

This is a non-hierarchical, three-grade system. 

The inital stock vector is: 



n_(0) = (300, 200, 100) 

The transition matrix is: 

( . b 0 .15 .05 

.75 .20 0 

.05 C .90 

and 

( .74 .164 .096 

.785 .2035 .0115 

.085 .0035 .9115 

The 0** vector is: 

0 * * = ( 1 , 0 , 0 ,) . 

The recruitment proportion vector is: 




r = (. 70, -07, .23) . 



Then 



I - Q = 



(I - Q) * = 



. 26 


-. 164 


-.096 


. 785 


. 797 


-.012 


. 085 


-. 004 


. 089 


with 


ones. 


gives 


1 


. 164 


-.096 \ 


1 


.797 


-.012 I 


1 


.004 


.089 / 
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and the inverse is 



( (I - Q) *)-i 



( .430 
-.611 
-4.888 



.091 .479 

1.127 -.516 

- . 9 8*1 5.869 



Using Eguation 2.9 



( .430 

-.611 

-4. 888 



.091 .479 

1.127 -.516 

-.981 5.869 



= ( . 43. .09, .48) 

is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, 100, 

350, and 400 can be found in Appendix D, confirming that in 
steady state there will be 4 35?, 95i, and 48% respectively, in 
categories 1,2, and 3. As the printout shows, steady state 
distribution (to two decimal accuracy) is achieved by 
t = 400 . 



2. Option 5: Multiplicative System Size 

This section addresses the option of multiplicative 
increase (decrease) in total system size for eacn time 
period. Starting with 

N (t) = 9N (t- 1 ) , for t = 1,2,... 
it follows that 

N (t) = 6*^(0), for t = 1,2,... 
from which M (t) can be determined as 

M(t) = N (t) — N (t— 1 ) = (6- 1) 0 r - IN (0) - 
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Then, letting (6- 1) 9 r_l N (0) replace M(t) in Equation 1.6, 
results in 

n (t) = n(t-1) Q +[ (0-1 N(0) ] r . (egn 2. 10) 

First, 6 > 0 must be true otherwise negative or zero system 
size may result. This leaves two cases to be examined. 

Case 1: 0 < 6 < 1 

In this case 

limi t N ( t) =0 

t — ► oo 



because as t — ► oo , 0*- i — ► 0 for 0 < 0 <1 . Finally, if 
0=1 then 0-1 = 0 and N (t) =0 for all t. Therefore, 

n. = limit n_ (t) = 0 

t oo 



must also hold. This means that such a system becomes empty 
in steady state. 

Case 2: 0 > 1 



Equation 2.10 can be rewritten as 
n (t) 



N (t) 



n ( t-1 ) (0-1) 9^-1 N(0) 

Q + r 



N(t) 



N{t) 



Since N (t) can be written as 0N(t-1) or as 0^(0) , the above 
equation gives 

n (t) n ( t- 1 ) 1 (0-1) e + ~iN (0) 

“ - - *Q ♦ x r 



N (t) 



N (t-1 ) 0 



G*N (0) 



or 



n (t) 



n ( t- 1 ) 1 

• - *Q + 



N (t) N (t-1 ) 0 

Taking the limit as t — *• «> results in 



e-i 

© 
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or 



where 



n = limit 
~ t - ~ 



~ 1 < 0 - 1 ) 

n = n - Q + r 

- ~ e 9 “ 



n [©I-Q] = ( 0— 1)_r 

n (t) 



N (t) 



is the steady state distribution vector. 

Since 6 > 1 in this case, the matrix 9I-Q has a 
unique inverse according to Theorem 2, so the result 



n = (9-1) r[ 91 -Q ] — 1 



( egn 2. 1 1) 



is obtained for the steady state distribution vector. After 
steady state has been reached it follows that the formula 

n (t) = 9 t N (0) n 

can be used tc find the stocks at time t. Two examples are 
worked out below to illustrate Equation 2.11 . 

E xample (a) 

This is a hierarchical, four category system. 

The inital stock vector is: 



P = 



(129, 


74, 28 , 


11) 


• 


is: 


/ .728 


. 102 


0 


° \ 


( 0 


.83 


.046 


0 I 


\ 0 


0 


.867 


. 033 1 


\ 0 


0 


0 


. 902 / 



<42 



and 



Q = 




.102 

,83 

0 

0 




The multiplicative factor is: 

9 = I.Ol’ . 

The recruitment proportion vector is: 



Then 



ei - q = 



and the inverse is 



(91 - Q)-» = 



0, 0, 0) 


* 






/ .112 


102 


0 


° \ 


/ -.124 


. 18 


-.046 


0 1 


l- 1 


0 


.143 


-.033 1 


\ -. 098 


0 


0 


. 108 / 



55.4 1 


'31.40 


10. 10 


3. 09 


51.04 


34.48 


11. 09 


3. 39 


50.35 


28.53 


16. 17 


4. 94 


50. 28 


28.49 


9. 17 


12. 06 



Using Eguation 2.11 



n = (-0 1 ) (1 ,0,0,0) 



55. 4 1 


3 1.40 


10. 10 


3.09 


51.04 


34.48 


1 1. 09 


3. 39 


50.35 


28.53 


16. 17 


4. 94 


50. 28 


28. 49 


9. 17 


1 2. 06 



= (. 5 54 r .314, .101, .031) 

is the steady state distribution vector. 
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The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, and 100, 
can be found in appendix E, confirming that in steady state 
there will be 55$, 31%, 10%, and 3% respectively, in catego- 
ries 1,2,3, and 4. As the printout shows, steady state 
distribution (to two decimal accuracy) is achieved by 
t = 50 . 

Example (b) 

This is a non-hierarchical, three-grade system. 



The inital stock vector is: 

n (0) = (300, 200, 1 00) 

The transition matrix is: 



and 



60 


. 15 


.05 \ 


75 


- 20 


° 


05 


0 


.90 / , 


.74 




164 .096 



Q = 



785 .2035 .0115 

085 .0035 .9115 



The multiplicative factor is: 

e = i.oi 

The recruitment proportion vector is: 

r = (. 70, .07 , . 23) . 



Then 



ei - Q 



/ . 27 

I -. 785 
\ -. 085 



-. 164 
. 8065 
-. 0035 



-. 096 
-.0115 

. 0985 
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and the inverse is 



(61 - Q)~ l 



/45. 32 9. 41 45. 27 

I 44.69 10.52 44.79 

\ 40. 70 8.50 50. 81 



Osing Equation 2.11 



n 



(.0 1 ) 



(-7,. 07, 



/ 4 5. 3 2 
( 44.69 
\ 40.70 



9.41 45.27 \ 

10.52 44.79 ] 

8. 50 50.81 / 



= (.44, . 09, .47) 

is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, and 100, 
can be found in Appendix E, confirming that in steady state 
there will oe 44??, 9 % , and 47*? respectively, ir. categories 
1,2, and 3. As the printout shows, steady state distribu- 
tion (tc two decimal accuracy) is achieved by t - 100 . 
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III. SOMMABY 



This chapter is a synopsis of the options that describe 
the five different conditions of growth. In all options, 
the restriction of w t - > 0 for all i = 1,2,. ..,k was assumed, 
since attrition must be allowed from any category in a 
manpower system. Options 1 through 3 employed the basic 
prediction eguation (Equation 1.2) and used constant, addi- 
tively and multiplicatively increasing (decreasing) recruit- 
ment, R(t), to derive the steady state behavior of a 

manpower system. For options 4 and 5, Eguation 1.6 used 
additively and multiplicatively increasing (decreasing) 
total system size, M (t) , to derive the steady state behavior 
of a manpower system. 

Derivation of the steady state stock vector or distribu- 
tion vector required some standard results of matrix algebra 
summarized in Theorems 1 and 2 of Chapter II. In addition 
to the analytical derivation of the steady state equations, 
examples have been worked out in each case (excepting some 
trivial cases) to illustrate the use of the analytical 
results. The formulas have been programmed in the AFL 
language and the program listing is given in Appendix F. 
This program (called a function in APL) has been integrated 
into the main APL program at NPS which forecasts future 
stocks of a manpower system using the Markov Chain theory 
summarized m Chapter I. Appendices A through E show the 
printouts of the terminal sessions where the examples 
referred to above are worked out to verify the steady state 
results obtained in each analytical example. 
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RPP ENDIX & 

COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 1 



1. Example (a) 



START 

DO YOU WISH TC ENTER DATA? 

0 NC 

1 Y ES 

« 

1 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

1 

ENTER N (INITIAL STOCK VECTOR) 

1 29 74 28 11 

ENTER THE PROMOTION RATE VECTOR. 

THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE FIRST 
3 CLASSES. 

. 102 . 046 .033 

ENTER THE WASTAGE RATE VECTOR. 

THE 4TH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAI 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

.17 .124 .1 .098 

ENTEP THE NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

1 

ENTER R (RECRUITMENT VECTOR) 

35 0 0 0 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

7 QUIT PROGRAM 

1 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 
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0 0 
0. 046 0 
0. 867 0.033 
0 0. 902 



P MATRIX 
0. 728 0. 102 
0 0.83 

0 0 

0 0 

N VECTOR 
129 74 28 1 1 
OPTION =1 
R VECTOR 
35 0 0 0 

WOULD YOU LIKE TO CHANGE ANY OF THE DATA? 

0 NO 

1 Y ES 

7 QUIT PRCGEAM 



0 

ENTER THE NUMBER OF THE YEAR 
1 0 



YOU WISH TO SEE 



DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

m 

o 



T 




N 


PERCENT R 


0 


1 


129 i 


( 53) 




2 


74 i 


31 




3 


28 | 


12) 




4 


1 1 1 


5) 




TOTAL 


242 i 


[ 100) 


10 


1 


129 | 


; 53) 




2 


77 i 


32 




3 


27 i 


11 




4 


10 | 


4) 




TOTAL 


242 i 


1 100) 



35 



DC YOU KISH TO SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



50 



0 



1 


129 | 


[ 53) 


2 


77 | 


32 


3 


27 | 


11) 


4 


9 i 


4) 


TOTAL 


242 | 


100) 



35 
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DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



100 



DO YOU 

0 NO 

1 Y ES 
7 QUIT 



1 


129 ( 


53' 


2 


77 \ 


32 ( 


3 


27 { 


1 1 


4 


9 


4’ 


TOTAL 


242 ( 


100 


WISH TO 


SEE ANY OTHER 


YEARS 



35 



PRCGEAM 



0 



DO YOU WISH TO 


SEE THE STEADY STATE 


VECTOR? 


0 NO 








1 Y ES 

• 








1 

999 1 


129 i 


( 5 3) 




2 


77 i 


32) 




3 


27 | 


11) 




4 


9 I 


4) 




TOTAL 


242 | 


100) 


35 



ARE YOU THROUGH? 

0 NC 

1 Y ES 

m 

i 



2- Exa m ple (b) 



START 

DC YOU WISH TC ENIER DATA? 

0 NC 

1 YES 

1 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

3 
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ENTER N (INITIAL STOCK VECTOR) 

300 200 100 

ENTER P (TRANSITION MATRIX) BY ROWS 
ENTER 1 TH ROW 

.6 .15 .05 

ENTER 2 TH ROW 

.75 .2 0 
ENTER 3 TH ROW 

. 05 0 .9 

ENTER THE NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

1 

ENTER R (RECRUITMENT VECTOR) 

70 7 23 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

7 QUIT PROGRAM 

1 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 
0. 6 0 . 1 5 0 . 05 

0. 75 0 .2 0 

0. 05 0 0.9 

N VECTOR 
300 200 100 
OPTION =1 
R VECTOR 
70 7 23 

WOULD YOU LIKE TO CHANGE ANY CF THE DATA? 

0 NO 

1 YES 

7 QUIT PRCGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 

1 0 

DC YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO - 

1 YES 

0 



50 



T N PERCENT E 



0 


1 


300 | 


[ 50) 




2 


200 < 


33) 




3 


100 | 


17) 




TOTAL 


600 | 


[ 100) 


10 


1 


351 | 


[ 48) 




2 


75 | 


10) 




3 


302 | 


41 




TOTAL 


728 < 


121 



DO YOU KISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



50 1 


375 ( 


43) 




2 


79 ( 


9) 




3 


415 ( 


48) 




TOTAL 


868 ( 


145) 


100 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

m 

i 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



100 1 


376 ( 


43) 


2 


75 \ 


9 


3 


418 ( 


48) 


TOTAL 


873 ( 


146 



DC YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 
0 

DO YOU WISH TC SEE THE STEADY STATE VECTOR? 

0 NC 

1 Y ES 

1 



51 



999 1 


376 < 


43) 




2 


79 | 


9) 




3 


418 i 


48 




TOTAL 


873 ( 


100) 


100 



ABE YOU THROUGH? 

0 NO 

1 YES 

« 

m 

i 
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APPENDIX B 



COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 2 



1 . Exa m ple (a) 



START 

DO YOU WISH TO ENTER DATA? 

0 NC 

1 Y ES 

1 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

1 

ENTER N (INITIAL STOCK VECTOR) 

129 74 28 11 

ENTER THE PROMOTION RATE VECTOR. 

THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE FIRST 
3 CLASSES. 

. 102 - 046 .033 

ENTER THE WASTAGE RATE VECTOR. 

THE 4TH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

.17 .124 .1 .098 

ENTER THE NUMBER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

2 

ENTER R (RECRUITMENT VECTOR) 

35 0 0 0 

ENTER ADDITIVE INCREASE 
5 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 

1 
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WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 

0. 728 0. 102 0 0 

0 0.83 0. 046 0 

0 0 0. 867 0.033 

0 0 0 0.902 

N VECTOR 
129 74 28 1 1 
OPTION =2 
INC =5 

TOTAL RECRUITMENT 
35 

RECRUITMENT PROPORTION VECTOR 
10 0 0 

WOULD YCU LIKE TO CHANGE ANY OF THE DATA? 

0 NO 

1 YES 

7 QUIT PROGRAM 

# 

0 

ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
1 0 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 



1 YES 
0 
T 


N 


PERCENT 


R 


0 1 


129 | 


; 53j 




2 


74 < 


31) 




3 


28 | 


12) 




4 


1 1 1 


5 




TOTAL 


242 | 


100) 




10 1 


248 | 


; 63 ) 




2 


106 { 


27 ) 




3 


29 | 


7 ) 




4 


10 


3 




TOTAL 


392 | 


162) 





DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DO YCU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



54 



280 



50 



1 


980 | 


58 


2 


523 j 


31 


3 


152 j 


9 


4 


39 i 


2 


TOTAL 


1695 < 


700 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 

# 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



100 



1 


1899 < 


' 56 


2 


1075 \ 


31 


3 


343 i 


10 


4 


102 < 


3 


TOTAL 


3419 { 


1413 
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DO YOU WISH TC SEE ANY OTHER YEARS? 

0 N C 

1 Y ES 

7 QUIT PROGRAM 
1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
9 00 

DO YOU WISH TO SEE THE INTERVENING YEAPS? 

0 NO 

1 YES 



0 

900 



1 

2 

3 

4 

TOTAL 



16605 | 


[ 54 


9898 | 


32 


3395 1 


1 1 


1130 < 


4 


31028 | 


[ 12821 



4530 



DC YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
9 75 

DO YOU WISH IC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 
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975 1 17984 < 53) 

2 10725 { 32) 

3 3681 ( 11) 

4 1226 ( 4) 

IOTA L 33617 ( 13891) 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

DO YOU WISH TC SEE THE STEADY STATE 

0 NC 

1 Y ES 

1 



PERCENTAGES ARE 


IN STEADY 


STATE 


999 1 


17984 ( 


53) 


2 


10725 


32) 


3 


3681 ( 


11) 


4 


1226 { 


4) 


TOTAL 


33617 ( 


100) 



ARE YOU THROUGH? 

0 NO 

1 YES 

1 



2. Exa mp le <b) 



START 

DO YOU WISH TC ENTER DATA? 

0 NO 

1 YES 

1 

ENTER THE NUMBER OR THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

3 

ENTEP N (INITIAL STOCK VECTOR) 

# 

300 200 100 

ENTER P (TRANSITION MATRIX) BY ROWS 
ENTER 1 TH ROW 

- 6 .15 .05 

ENTER 2 TH ROW 

.75 .2 0 
ENTER 3 TH ROW 
. 05 0 .9 



4905 



VECTOR? 



4905 
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ENTER THE NUMBER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

* 

2 

ENTER R (RECRUITMENT VECTOR) 

70 7 23 

ENTER ADDITIVE INCREASE 
25 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

7 QUIT PROGRAM 

* 

1 

NOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 

0.6 0.15 0.05 

0.75 0.2 0 

0. 05 0 0.9 

N VECTOR 

300 200 100 

OPTION =2 

INC =25 

TOTAL RECRUITMENT 

100 

RECRUITMENT PROPORTION VECTOR 
0.7 0.07 0. 23 

WOULD YOU LIKE TO CHANGE ANY OF THE DATA? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 0 
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DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

• 

0 

T N PERCENT R 



1 

2 

3 

TOTAL 


300 ( 
200 ( 
100 ( 
600 ( 


501 

33) 

17) 

100) 


1 


801 ( 


53) 


2 


162 ( 


11 ) 


3 


561 ( 


37 


TOTAL 


1523 


254) 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 



325 



1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DC YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



50 


1 


4374 ( 


46) 




2 


914 ( 


10 




3 


4260 ( 


45 




TOTAL 


9548 ( 


1591) 


YOU 


WISH TC 


SEE ANY OTHER 


YEARS? 



0 NO 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 

100 1 9066 ( 

2 1903 \ 

3 9469 < 

TOTAL 20438 ( 



44 

9 

46 

3406 



2575 



DO YOU WISH TC SEE ANY OTHER YEARS?' 

0 NO 

1 YES 

7 QUIT PROGRAM 



1 
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ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
3 50 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



350 



1 

2 

3 

TOTAL 



32550 ( 43) 

0853 ( 9) 

35585 ( 47) 

74987 { 12498) 



8825 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU FISH TO SEE 
4 00 



DO YOU WISH TC 

0 NO 

1 YES 


SEE THE INTERVENING YEARS? 


0 


400 1 


37246 | 


43) 


2 


7843 1 


9) 


3 


40808 


46) 


TOTAL 


85897 


14316) 10075 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

DO YOU WISH TC SEE THE STEADY STATE VECTOR? 

0 NC 

1 Y ES 

1 

PERCENTAGES ARE IN STEADY STATE 



1 


37246 ( 


43) 




2 


7843 ( 


9} 




3 


40808 ( 


48) 




TOTAL 


85897 ( 


100) 


10075 



ARE YOU THROUGH? 

0 NC 

1 Y ES 

1 
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APPENDI X C 

COMPUTER PRINTOUTS OF EXAMPLES FOE OPTION 3 



1 . E xample (a ) 



S1AET 

DO YOU WISH TC ENTER DATA? 

0 NC 

1 Y ES 

1 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

1 

ENTER N (INITIAL STOCK VECTOR) 

1 29 74 28 11 

ENTER THE PROMOTION RATE VECTOR. 

THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE FIRST 
3 CLASSES. 

. 102 . 046 .033 

ENTER THE WASTAGE RATE VECTOR. 

THE 4TH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAI 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

.17 .124 .1 .098 

ENTER the NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

3 

ENTER R (RECRUITMENT VECTOR) 

35 0 0 0 

ENTER MULTIPLICATIVE FACTOR 

1 .01 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 

1 
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WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 

0. 728 0.102 0 0 

0 0.83 0. 046 0 

0 0 0. 867 0.033 

0 0 0 0. 902 

N VECTOR 
129 74 28 1 1 
OPTION =3 
FACT =1.01 

TOTAL RECRUITMENT 
35 

RECRUITMENT PROPORTION VECTOR 
10 0 0 

WOULD YOU LIKE TO CHANGE ANY CF THE DATA? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 0 



DO YOU WISH TC 

0 NO 

1 YES 

0 

T 


SEE THE INTERVENING 
N PERCENT 


YEARS 
• R 


0 1 


129 I 


[ 53] 


1 


2 


74 I 


[ 31 




3 


28 i 


12; 




4 


1 1 1 


5' 




TOTAL 


242 | 


[ 100) 


10 1 


137 I 


54] 


1 


2 


79 | 


3 1 




3 


27 | 


11] 




4 


10 


; 4 


| 


TOTAL 


253 i 


105 





DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PRCG RAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 
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50 



1 


204 | 


[ 55 


2 


116 < 


31 


3 


37 | 


10 


4 


1 1 1 


3 


TOTAL 


368 | 


152 
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DC YOU WISH TC SEE ANY OTHER 

0 NO 

1 YES 

7 QUIT PROGRAM 



YEARS? 



1 

ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 

100 



94 



1 


336 ( 


55 


2 


190 \ 


31 


3 


61 


10 


4 


19 ( 


3 


TOTAL 


606 


250 



DC YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PRCGEAM 



0 

DO YOU WISH 

0 NC 

1 Y ES 

1 



TO SEE THE STEADY STATE VECTOR? 



PERCENTAGES ARE 


IN STEADY 


STATE 


999 1 


336 ( 


55 


2 


190 ( 


31 


3 


61 { 


10 


4 


19 ( 


3 


TOTAL 


606 ( 


100 
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ARE YOU THROUGH? 

0 NO 

1 YES 



1 



2- Exa mple <b) 

START 

DO YOU WISH TO ENTER DATA? 

0 NO 

1 YES 

1 
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ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

3 

ENTER N {INITIAL STOCK VECTOR) 

300 200 100 

ENTER P (TRANSITION MATRIX) BY ROWS 
ENTER 1 TH ROW 

.6.15 .05 
ENTER 2 TH ROW 

.75 .2 0 
ENTER 3 TH ROW 

. 05 0 .9 

ENTER THE NUMBER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

3 

ENTER R (RECRUITMENT VECTOR) 

70 7 23 

ENTER MULTIPLICATIVE FACTOR 

1 .01 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

7 QUIT PROGRAM 

1 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 



0.6 0.15 0. 05 

0. 75 0.2 0 

0.05 0 0.9 

N VECTOR 
300 200 100 

OPTION =3 

FACT =1.0 1 

TOTAL RECRUITMENT 

100 

RECRUITMENT PROPORTION VECTOR 
0.7 0.07 0. 23 
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WOULD YOU LIKE 10 CHANGE ANY OF THE DATA? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 

1 0 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

o 

T N PERCENT R 



0 1 


300 ( 


50) 


2 


200 ( 


33 


3 


100 ( 


17 


TOTAL 


600 { 


100) 


10 1 


370 ( 


49) 


2 


78 


10 


3 


312 


41 


TOTAL 


761 ( 


127) 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU KISH TO SEE 
50 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



50 1 


575 ( 


44 


2 


121 { 


• 9 


3 


603 ( 


46 


TOTAL 


1298 < 


216 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 
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100 1 


947 | 


' 44) 


2 


199 1 


9 


3 


996 | 


47 


TOTAL 


2141 I 


357 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

DO YOU WISH TO SEE THE STEADY STATE VECTOR? 

0 NO 

1 Y ES 

1 

PERCENTAGES ARE IN STEADY STATE 



999 1 


947 < 


44) 




2 


199 ( 


9) 




3 


996 ( 


47) 




TOTAL 


214 1 


100) 


268 



ARE YOU THROUGH? 

0 NO 

1 Y ES 

1 
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COMPUTER 



APPENDIX D 

PRINTOUTS OF EXAMPLES FOR OPTION 4 



1. Case J Exam ple (a) 



START 

DO YOU WISH TC ENTER DATA? 

0 NO 

1 YES 

1 

ENTER THE NUMEER OF THE MODEL IYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

1 

ENTER N (INITIAL STOCK VECTOR) 

1 29 74 28 11 

ENTER THE PROMOTION RATE VECTOR. 

THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE FIRST 
3 CLASSES. 

. 102 . 046 .033 

ENTER THE WASTAGE RATE VECTOR. 

THE 4TH VALUE SHOULD INCLUDE THE TOTAL RATE FOE THAT 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

.17 .124 .1 .098 

ENTER THE NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

4 

ENTER R PROP (RECRUITMENT PROPORTION VECTOR) 

10 0 0 

ENTER ADDITIVE INCREASE 
0 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

7 QUIT PROGRAM 

1 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 



66 



P MATRIX 

0. 728 0. 102 0 0 

0 0.83 0. C46 0 

0 0 0. 867 0.033 

0 0 0 0. 902 

N VECTOR 
129 74 28 1 1 
OPTION =4 
INC =0 

RECRUITMENT PROPORTION VECTOR 
10 0 0 

WOULD YOU LIKE TO CHANGE ANY CF THE DATA? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 0 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



T 


N 


PERCENT R 


0 1 


129 | 


[ 53) 


2 


74 | 


31 


3 


28 1 


12) 


4 


1 1 i 


5 


TOTAL 


242 | 


[ 100) 


10 1 


129 I 


f 53) 


2 


77 \ 


32) 


3 


27 | 


1 1) 


4 


10 i 


4) 


TOTAL 


242 i 


100 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DC YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 XES 

0 



1 


125 | 


' 53) 


2 


77 i 


32) 


3 


27 ( 


11) 


4 


9 i 


4) 


TOTAL 


242 | 


100) 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 



1 

ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



0 

1 

7 



100 


1 


129 ( 


53 




2 


77 { 


32 




3 


27 ( 


11 




4 


9 ( 


4 




TOTAL 


242 ( 


100 


YOU 


WISH TC 


SEE ANY OTHER 


YEARS 



35 



NC 
Y ES 

QUIT PROGRAM 



0 



DO YOU WISH 

0 NC 

1 Y ES 

1 

PERCENTAGES 
999 



TO SEE THE STEADY STATE VECTOR? 



ARE IN STEADY STATE 



1 


129 i 


f 53) 


2 


77 { 


32) 


3 


27 | 


11) 


4 


9 { 


4) 


TOTAL 


242 i 


100 



35 



ARE YOU THROUGH? 

0 NC 

1 YES 

1 



2. Case J Example ( b) 



START 

DO YOU WISH TO ENTER DATA? 

0 NC 

1 YES 

1 • 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

3 
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ENTER N (INITIAL STOCK VECTOR) 

300 200 100 

ENTER P (TRANSITION MATRIX) EY RONS 
ENTER 1 TH ROW 

.6.15 .05 
ENTER 2 TH ROW 

.75 -.2 0 
ENTER 3 TH ROW 

. 05 0 .9 

ENTER THE NUMEEE OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

4 

ENTER R PROP (RECRUITMENT PROPORTION VECTOR) 

.7.07 .23 

ENTER ADDITIVE INCREASE 
0 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

7 QUIT PROGRAM 

1 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 
0. 6 0.15 0.C5 

0. 75 0.2 0 

0. 05 0 0.9 

N VECTOR 
300 200 100 
OPTION =4 
INC =0 

RECRUITMENT PROPORTION VECTOR 
0.7 0.07 0. 23 



WOULD YOU LIKE TO CHANGE ANY OF THE DATA? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 0 
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DO you WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



T 




N 


PERCENT 


0 


1 


300 | 


50) 




2 


200 < 


33'i 




3 


100 < 


17 




TOTAL 


600 < 


: ioo) 


10 


1 


284 i 


f 47) 




2 


61 | 


10) 




3 


255 i 


43) 




TOTAL 


600 l 


1 00 i 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 
1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



50 1 


258 j 


[ 43) 




2 


54 i 


9) 




3 


287 { 


48) 




TOTAL 


600 i 


100) 
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DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

« 

1 

ENTER THE NUMEER OF THE YEAR YOU KISH TO SEE 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



100 1 


258 | 


4 3 


2 


54 i 


9 


3 


287 | 


48 


TOTAL 


600 | 


100 
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DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 
0 

DO YOU WISH TO SEE THE STEADY STATE VECTOR? 

0 NC 

1 Y ES 

m 

i 

PERCENTAGES ARE IN STEADY STATE 



999 1 


258 ( 


43) 


2 


54 ( 


9' 


3 


287 | 


48) 


TOTAL 


600 | 


100) 



ARE YOU THROUGH? 

0 NC 

1 Y ES 

1 



3. Case 2 Examp l e (a) 



START 

DO YOU WISH IC ENTER DATA? 

0 NC 

1 YES 

1 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

1 

ENTER N (INITIAL STOCK VECTOR) 

1 29 74 28 11 



ENTER THE PROMOTION RATE VECTOR. 

THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE FIRST 
3 CLASSES. 

. 102 . 046 .033 

ENTER THE WASTAGE RATE VECTOR. 

THE 4TH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

.17 .124 .1 .098 

ENTER THE NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

4 
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ENTER R PROP (RECRUITMENT PROPORTION VECTOR) 

10 0 0 

ENTER ADDITIVE INCREASE 
5 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 

m 

i 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 

0. 728 0. 102 0 0 

0 0.83 0. 046 0 

0 0 0. 867 0.033 

0 0 0 0.902 

N VECTOR 
129 74 28 1 1 

OPTION =4 

INC =5 

RECRUITMENT PROPORTION VECTOR 
10 0 0 

WOULD YOU LIKE TC CHANGE ANY OF THE DATA? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 

1 0 



DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 

T N PERCENT R 



1 


129 I 


’ 53) 


2 


74 i 


31) 


3 


28 1 


12 


4 


1 1 i 


5 


TOTAL 


242 | 


[ 100) . 


1 


165 | 


( 57) 


2 


89 | 


30 


3 


28 ( 


10) 


4 


10 i 


[ 3) 


TOTAL 


292 i 


121) 
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DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 

DO YOO WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



50 



0 



1 


273 | 


: 56) 


2 


155 | 


31 


3 


49 | 


10) 


4 


15 ( 


3) 


TOTAL 


492 


203) 
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DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



100 



0 



1 

2 

3 

4 

TOTAL 



407 

235 

77 

24 

742 



III 

307) 



113 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

«• 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
9 00 



DO YOU WISH TC 

0 NO * 

1 YES 


SEE THE INTERVENING 


0 


900 1 


2537 | 


; 54) 


2 


1513 i 


32) 


3 


519 l 


1 1 


4 


173 i 


4) 


TOTAL 


4742 | 


1960) 



692 
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DO YOU HISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 

i 

ENTER THE NUMEER OF THE YEAR YOU HISH TO SEE 
9 50 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



53) 

32) 

1 1 ) 

4) 

2063) 728 



950 



1 

2 

3 

4 

TOTAL 



2670 

1593 

547 

182 

4992 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 
0 

DO YOU WISH TC SEE THE STEADY STATE VECTOR? 

0 NC 

1 Y ES 

1 



percentages ape 


IN STEADY 


STATE 


959 1 


2670 ( 


53 


2 


1593 ( 


32 


3 


547 ( 


1 1 


4 


182 ( 


4 


TOTAL 


4992 ( 


100 



ARE YOU THROUGH? 

0 NO 

1 YES 

1 



4. Case 2 Exam pl e (b) 



START 



DO YOU WISH TC ENTER DATA? 

0 NC 

1 YES 

1 



74 



ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

3 

ENTER N (INITIAL STOCK VECTOR) 

300 200 100 

ENTER P (TRANSITION MATRIX) BY ROWS 
ENTER 1 TH ROW 

.6 .15 .0 5 

ENTER 2 TH ROW 

.75 .2 0 
ENTER 3 TH ROW 

. 05 0 .9 

ENTER THE NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

4 

ENTER R PROP (RECRUITMENT PROPORTION VECTOR) 
.7 .07 .23 

ENTER ADDITIVE INCREASE 
25 



ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 

m 

i 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 
0. 6 0 . 1 5 0 . 05 

0. 75 0.2 0 

0. 05 0 0.9 

N VECTOR 
300 200 100 

OPTION =4 
INC =25 

RECRUITMENT PROPORTION VECTOR 
0.7 0.07 0. 23 
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WOULD YOU LIKE 10 CHANGE ANY OF THE DATA? 

0 NO 

1 YES 

7 QUIT PROGRAM 
0 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 0 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 



YES 

0 

T 


N 


PERCENT 


R 


0 1 


300 ( 


' 50] 


1 


2 


200 | 


33 




3 


100 | 


17 




IOTA L 


600 1 


[ 100 1 




10 1 


417 | 


[ 491 




2 


88 1 


10 




3 


345 i 


41 




TOTAL 


850 | 


142 





DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 
1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 



DO YOU WISH IC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



50 1 


827 | 


45) 


2 


173 


9 


3 


849 | 


46 


TOTAL 


1850 < 


303 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 
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100 1 


1365 ( 


44) 




2 


287 ( 


s l 




3 


1448 ( 


47 




TOTAL 


3100 ( 


517 


382 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 



1 



ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
3 50 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



350 



1 

2 

3 

TOTAL 



4056 


f 43) 




854 


1 9) 




4440 


( 47) 




9350 


( 1558) 


1098 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 YES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
4 00 



DC YOU WISH TO 

0 NO 

1 YES 

0 

400 1 

2 
3 

TOTAL 



SEE THE INTERVENING YEARS? 



4594 ( 
967 / 
5039 ( 
IOoOO { 



48) 

1767 



1241 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PR CG RAM 
0 

DO YOU WISH TO SEE THE STEADY STATE VECTOR? 

0 NO 

1 Y ES 

1 
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PERCENTAGES ARE IN STEADY ST A IE 



999 1 


4594 ( 


43) 


2 


967 < 


9 


3 


5039 | 


48) 


TOTAL 


10600 { 


100) 



ARE YOU THROUGH? 

0 NO 

1 Y ES 

1 



1241 
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APP ENDI X E 

COMPUTER PRI8TOOTS OF EXAMPLES FOE OPTION 5 



1. Exa m ple (a) 



S1ART 

DO YOO WISH TO ENTER DATA? 

0 NC 

1 Y ES 

1 

ENTER THE NUMEER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

• 

1 

ENTER N (INITIAL STOCK VECIOR) 

1 29 7 4 28 11 

ENTER THE PROMOTION RATE VECTOR. 

THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE 
3 CLASSES. 

. 102 . 046 .033 

ENTER THE WASTAGE FATE VECTOR. 

THE 4TH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

17 12 4 1 0 98 

ENTER THE NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

5 

ENTER R PROP (RECRUITMENT PROPORTION VECTOR) 

10 0 0 

ENTER MULTIPLICATIVE FACTOR 

1 .01 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 

1 



FIRST 
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WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 



1 



P MATRIX 
0. 728 0. 102 
0 0.83 

0 0 
0 0 
N VECTOR 
129 74 28 1 1 
OPTION =5 
FACT =1.0 1 
RECRUITMENT 
10 0 0 



0 

0. 046 
0. 867 
0 



0 

0 

0.033 

0.902 



PROPORTION VECTOR 



WOULD YOU LIKE TC CHANGE ANY OF THE DATA? 

0 NC 

1 YES 

7 QUIT PROGRAM 



0 



ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 

1 0 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 

T N PERCENT E 



1 


129 ( 


: 53) 


2 


74 | 


31 


3 


28 | 


12) 


4 


1 1 1 


[ 5 


TOTAL 


242 ( 


: i oo) 


1 


148 | 


[ 55) 


2 


82 j 


31 


3 


27 < 


10) 


4 


10 i 


1 V 


TOTAL 


268 1 


111) 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 

1 

ENTER THE NUMEER OF THE YEAR YOU FISH TO SEE 
50 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



90 



1 


222 ( 


55 


2 


125 ( 


31 


3 


40 ( 


10 


4 


12 ( 


3 


TOTAL 


400 


165 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 



1 

ENTER THE NUMEER OF THE YEAR YOU WISH 
1 00 

DO YOU WISH TC SEE THE INTERVENING YEA 

0 NO 

1 YES 

0 



100 



1 

2 

3 

4 

TOTAL 



362 ( 
206 
66 ( 
20 { 
654 



ifl 

270) 



TO SEE 
RS? 



101 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 v p s: 

7 QUIT PROGRAM 

m 

n 



DO YOU WISH TO SEE THE STEADY STATE VECTOR? 
C NC 
1 Y ES 

1 



PERCENT 

999 



AGES ARE IN STEADY STATE 



1 


362 i 


( 55) 


n 


206 ( 


: 31) 


3 


o6 1 


[ 10 


4 


20 | 


3) 


TOTAL 


654 | 


100 



ARE YOU THROUGH? 

0 NO 

1 YES 

1 



10 1 



2. Example (b) 

START 

DO YOU WISH TC ENTER DATA? 

0 NC 

1 YES 

1 
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ENTER THE N OH EER OF THE MODEL TYPE 

1 MARKOV HIERARCHICAL 

2 MARKOV LENGTH OF SERVICE 

3 MARKOV GENERAL 

4 VACANCY 

3 

ENTER N (INITIAL STOCK VECTOR) 

3 00 2 00 1 00 

ENTER ? (TRANSITION MATRIX) BY ROWS 
ENTER 1 TH ROW 

.6.15 .05 

ENTER 2 TH ROW 

.75 .2 0 
ENTER 3 TH ROW 

05 0 .9 

ENTER the NUMEER OF THE RECRUIT TYPE 

1 FIXED RECRUIT VECTOR 

2 ADDITIVE (RECRUIT SIZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

u ADDITIVE (SYSTEM 31 Z*) 

5 1UITI PLICaTI VE (SYSTEM SIZE* 

5 

ENTER R PROP (RECRUITMENT PROPORTION VECTOR) 

.7 .07 .23 

ENTER MULTIPLICATIVE FACTOR 

1 .01 

ENTER THE PERCENT CODE 

0 NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADS SIZE 

7 QUIT PROGRAM 

1 

WOULD YOU LIKE TO SEE THE ENTERED DATA? 

0 NO 

1 YES 

1 

P MATRIX 
0.6 0.15 0.05 

0. 75 0 .2 0 

0. 05 0 0.9 

N VECTOR 
300 200 100 
OPTION =5 
FACT =1.01 

RECRUITMENT PROPORTION VECTOR 
0.7 0. 07 0. 23 

WOULD YOU LIKE TO CHANGE ANY OF THE- DATA? 

0 NC 

1 Y ES 

7 QUIT PROGRAM 
0 
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ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 

1 0 

DO YOU WISH TC SEE THE INTERVENING YEARS? 

0 NO 

1 YES 



0 



T 




N 


PERCENT R 


0 


1 


300 | 


; 50) 




2 


200 l 


33 




3 


100 i 


17) 




TOTAL 


600 i 


[ 100) 


10 


1 


316 i 


48) 




2 


67 | 


10 




3 


277 i 


42) 




TOTAL 


66 1 | 


1 10) 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NC 

1 Y ES 

7 QUIT PROGRAM ' 

1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
50 



DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



1 


435 ( 


44 


9 


91 \ 


9 


3 


458 ( 


46 


TOTAL 


985 ( 


164 



DO YOU WISH TC SEE ANY OTHER YEARS? 

0 NC 

1 YES 

7 QUIT PROGRAM 



123 



1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
1 00 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

0 NO 

1 YES 

0 



1 


717 ( 


44 


2 


150 \ 


9 


3 


754 


47 


TOTAL 


1621 


270 



83 



DO YOU WISH TO SEE ANY OTHER YEARS? 

0 NO 

1 Y ES 

7 QUIT PROGRAM 
0 

DO YOU WISH TC SEE THE STEADY STATE VECTOR? 

0 NC 

1 Y ES 

1 

PERCENTAGES ARE IN STEADY STATE 



999 1 


7 17 { 


44) 




2 


150 ( 


9 




3 


754 ( 


47) 




TOTAL 


1621 ( 


100) 


203 



ARE YOU THROUGH? 

0 NC 

1 YES 

1 
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APP ENDIX P 

APL COMPUTER PROGRAM FOR STEADY STATE DISTRIBUTION FUNCTION 

This is the APL function that is used to find the steady 
state distribution vector for all options. 



1 

2 

3 

4 



7 ' 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 



V STDST ; Q ; I DENT ;RV-,W; NSS ; NHAT 1 ; NEAT 2 ; FACTI ; IQ ; ZERO . 
N*-N J 

RV*- (1,K) pRPROP 
I DENT*- (K,K)p(l,KpO) 

W*-(l,K)ol- ( +/ ( 2 ] P) 

Q*-P+ ( (W) + .*RV) 

( TYPE= 1 2 3 4 5 ) / FIX , ADDREC .MULTREC , ADDSYS , MULTSYS 
FIX : N*-R+ . xg (I DENT- P ) 

NSS*-N+ (+/N) 

ROUND - N 

+END 

ADDREC : NHAT1+RPROP+ . xg (IDENT-P) 

NSS+NHATH ( +/NHAT1 ) 

+END 

MULTREC : FACTI*-FACT*IDENT 
NHAT2*-R+ . xg ( FACTI- P ) 

NSS+NHAT2±(+/NHAT2 ) 

+END 

ADDSYS : IQ*-IDENT- Q 
IQ l ; 1 ] *-Kp 1 
ZERO *- 1 , (K- 1 loO 
NSS+ZERO+.x&IQ 
^END 

MULTSYS : FA CTI *-FA CT* I DENT 
NSS*- (FACT- 1 ) *RPR0P+ . *1 ( FACTI- Q ) 

END~ TCOUNT*-§99 
PERCENTS (0 . 5+100 *NSS) 

TOTPERCENT *- 10 0 
OUTPUT 
+ 0 
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